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Abstract 

In this paper, we propose a distributed rate allocation algorithm that minimizes the average decoding delay in inter-session 
network coding systems. We consider a scenario where network nodes, or peers, are organized in a mesh network, with every 
peer interested in the content of one of the several available sources. Our novel algorithm allows peers to determine the coding 
operations and the packet rates to be requested from the parent peers. A rate allocation problem is solved by every peer, which 
seeks the rates that minimize the average decoding delay for its children and for itself. To solve this non-convex optimization 
problem, we introduce the concept of equivalent packet flows, which permits to estimate the expected number of packets that every 
peer needs to collect for decoding. We then decompose our original rate allocation problem into a set of convex subproblems, 
which are eventually combined to obtain an approximate solution to the delay minimization problem. The results demonstrate 
II \ that the proposed scheme eliminates effectively the bottlenecks and reduces the decoding delay experienced by users with limited 

- , bandwidth resources. 
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c/3 ■ I. Introduction 

^ Distributed network architectures and protocols have gained much popularity over the past few years due to their scalability 
properties. Deployed initially for file sharing, today distributed systems, such as peer-to-peer (P2P) networks, are exploited for 
more demanding network applications such as live streaming, VoD, multi-party conferencing etc. The essential advantage of 
these systems over the traditional client-server architecture is their ability to sustain a large number of users without increasing 
■ the server load, as peers contribute their upload bandwidth to the system. This, however, comes at the cost of dynamic and 
unpredictable behavior of the peers, which renders the centralized routing challenging and necessitates distributed algorithms 
'/"J for data retrieval. In this context, network coding [1] has been considered recently as a solution to improve the performance 
of distributed systems. It removes the need for content reconciliation among the peers and offers decentralized control as well 
T— I as efficient adaptation to bandwidth variations and losses. 

A broad spectrum of distributed algorithms that utilize network coding exists in the literature. These works mainly focus on 
[ the case where a single data source from one or more servers is delivered to multiple users. It is common, however, that the 
I> \ network resources need to be shared by simultaneous applications. In such settings, inter-session network coding ||2l arises as 
k> a natural extension of network coding for efficient use of network resources with multiple concurrent sessions. Yet, the design 
_ of the network codes is not a trivial task; random mixing of all the sessions that exist in the network may lead to significant 
d • increase in the delay for a user to decode its source of interest from the combinations of different sessions. 

In this paper, we build on our previous work [3 1 and address the problem of designing a distributed rate allocation algorithm 
that decides how many packets of each session combination should be transmitted on the network links. We consider a scenario 
with multiple concurrent sessions that transmit data to peers organized in mesh network. The proposed protocol is receiver- 
driven and comprises two steps. First, the peer requests and receives information about its local neighborhood that includes its 
parents and children peers. Second, the peer requests from its parents intra- and inter-session network coded packets at specific 
rates. These rates are obtained by solving an optimization problem that seeks for the optimal rate allocation among different 
packet combinations. The objective of the distributed optimization algorithm is to minimize the average delay necessary for 
the peer and its children peers to decode their data of interest. The inclusion of the children peers in the delay minimization 
aims at avoiding selfish behaviors and provides incentives for peers to trade off part of their resources for an improvement in 
the overall network performance. 

The delay minimization problem is a priori non-convex. We approximate it with a set of convex subproblems by introducing 
the new concept of equivalent flows. An equivalent flow is defined for every component session of an inter-session combination. 
It can be regarded as a hypothetical flow with rate equal to the receiving innovative rate for the particular component session. 
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This leads to an estimation of the expected number of packets necessary for decoding a source of interest from a particular packet 
combination. Based on the equivalent flows representation, the original optimization problem is decomposed into several convex 
rate allocation subproblems that are easily solvable. Their solutions are then combined to yield an approximate yet effective 
solution to the optimal rate allocation. Simulation results demonstrate that the proposed scheme ehminates the bottlenecks and 
reduces the decoding delay experienced by users with limited resources. 
In summary, the contributions in this paper mainly consist in: 

• the formulation of a decoding delay optimization problem for inter-session network coding in peer-to-peer configurations, 

• the introduction of the novel concept of equivalent flows for approximate delay computation in inter-session network 
coding scenarios, 

• the design of a new distributed rate allocation algorithm for minimizing the average decoding delay. 

The rest of the paper is organized as follows. In Section we discuss the related work. We describe the scenario that 
we consider and the communication protocol in Section |lll] In the same section, we formulate the distributed rate allocation 
problem with inter-session network coding. The concept of equivalent flows is introduced in Section|lV] Our proposed distributed 
algorithm for delay minimal rate allocation is presented in Section[Vl In Section|VT]we evaluate the performance of the proposed 
rate allocation scheme and Section IVIII concludes our work. 

II. Related work 

With the recent advances in network coding research, the potentials of network coding have developed in the framework of 
P2P and overlay data delivery networks For example, Wang et al. fSl have proposed a design called that combines the 
random push strategy with random network coding. is shown to improve the performance of live streaming in terms of initial 
buffering delays, resilience to peer dynamics and reduced bandwidth load on the servers. The work in ||6] provides an analysis 
of the rate-delay-reliability trade-offs in a P2P streaming system. Random linear network coding is used to avoid duplicate 
packet reception. The authors derive upper and lower bounds on the minimum initial buffering required so that the playback 
interruption probability remains below a certain level. Network coding has also been considered for unequal error protection 
in overlay streaming systems as in fTl, where the authors propose a distributed receiver-driven algorithm for prioritized media 
delivery. 

Apart from the single session streaming cases, network coding has been also considered for multiple concurrent unicast and 
multicast scenarios. It has been shown that linear network coding is not sufficient for achieving the capacity bound |8 1; however, 
significant throughput gains can still be obtained with linear inter-session network coding as shown in [9|, which describes 
an implementation of opportunistic network coding for multiple unicast flows over wireless networks. Recently, several inter- 
session network coding algorithms have been proposed, mainly for data delivery in wireless networks [10], [TT], [T3l, 
|[T4|. Some of the works extend the COPE architecture [9| by considering application-specific features when designing the 
network codes. The work in lITOl for example studies the benefits of delaying packets at intermediate nodes in order to create 
more network coding opportunities. The proposed network coding scheme builds on COPE and incorporates an optimization 
framework that seeks for the optimal code and transmission policies that optimize the rate-distortion function. The performance 
of COPE and COPE-based systems degrades significantly in the presence of losses and network coding is turned off when 
the packet loss rate reaches a certain threshold. To deal efficiently with the packet losses, the authors in |12| propose a joint 
application of intra-session and inter-session network coding. Intra-session network coding is used for protection against packet 
losses, whereas inter-session network coding increases the throughput of the network. In order to characterize the capacity 
achieved with inter-session network coding for the 2-hop relay networks in the presence of losses, a flow based analysis is 
presented in lfT3l . The key idea is to regard packets as members of flows and not as independent entities as in ||9l, ifTOI . lfT2ll . 
A different approach for finding the feasible rate region is built on virtual multicasts ifT?! . The flow -based problem formulation 
stated in |14| provides a rate region which is at least as large as the rate region that can be achieved without inter-session 
network coding. 

While the benefits of inter-session network coding are well understood in the wireless scenarios, in wireline networks the 
construction of practical inter-session network coding algorithms is more challenging. The reason lies in the difference between 
the two communication media. The broadcast nature of wireless channels promotes the application of inter-session network 
coding through overhearing JP], lITOl . i.e., packets that are required for decoding can be overheard without wasting additional 
resources and decoding can be performed at every hop. This is not the case in wireline networks. Various theoretical aspects of 
inter-session network coding, such as sufficiency of linear codes and complexity of identifying coding opportunities, are studied 
in ifTSlI for the special case of pairwise coding in wireline networks. Kim et al. ||T6| propose a more generic solution that utilizes 
linear network coding and does not restrict the codes to specific classes such as pairwise or XOR coding. The coding strategy 
is determined with the help of Genetic Algorithms that optimize a certain cost objective. The work in [17] provides a different 
perspective on the design of inter-session network coding algorithms by exploiting the queue-length information to make the 
scheduling-routing-coding decisions. In |18|, the authors propose a low-complexity receiver-driven P2P system for delivery 
of multiple description coded data. The proposed scheme combines Raptor codes [191 with intra- and inter-session network 
coding in order to efficiently exploit the path diversity in the streaming overlay. An interesting application of inter- session 
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Fig. 1. Illustration of a multi-session scenario. Each source provides different data to the network. The users are organized in a mesh network, where each 
user requests a specific source data. 

network coding is further presented in ll20l . The authors propose a new gesture broadcast protocol and discuss the effectiveness 
of using inter-session network coding for application with low bit rate bursty streams, such as streams of gestures. 

To the best of our knowledge, there is however no work in the literature that addresses the problem of minimizing the 
average decoding delay in wireline mesh networks by distributed rate allocation in inter-session network coding. 

III. Data delivery with inter-session network coding 

A. Framework 

We consider a set of sources S and a set of peer nodes J\f that request data from different sources. We consider that the 
source data is segmented into blocks of Ng, {s E S) packets, and that the sources transmit simultaneously at rate Us, {s G S). 
The peers are organized in a wireline mesh network. The network is assumed to be directed and free of cycles. Thus, it can 
be modeled as a directed acyclic graph Q = (V, £), where V represents the set of network nodes, such that V = S U J\f, and 
£ is the set of connecting links between the network nodes. The directed link connecting any two nodes i and j is denoted 
as G £■ It is characterized by the link capacity 6y expressed in packets/sec and the average packet loss probability tt^. 
Furthermore, we consider that the servers have limited upload capacity while the peer nodes have limits on both the upload 
and download bandwidth, denoted as , i G V, and Cf, i G J\f, respectively. Finally, if nodes i and j are connected with the 
directed link we call node j as a child of node i, and node i is called the parent of node j. An example of such mesh 

network is illustrated in Fig. [T] 

The peer nodes act both as end users and relay nodes. Each peer is interested in receiving only one of the source dataQ 
Since the upload bandwidth of the sources is limited and only a small number of peers can acquire the requested packets 
directly from the sources, the majority of the peers are served by their parent peers. This implies that a peer may request and 
forward not only packets of the source that it has subscribed to, but also packets that are useful for its children peers. In order 
to increase the network throughput and alleviate the bottlenecks created by the limited network resources, we propose to allow 
the peer nodes to implement inter-session network coding. Inter-session network coding ||2| is an extension of network coding 
HJ to the case of multiple concurrent sessions (data sources) that share the same network resources. It essentially consists 
in combining packets from different sessions (sources), contrarily to intra-session network coding where only packets of the 
same session (source) participate in the packet combinations. When linear operations are considered, an inter-session network 
coded packet can be formally represented as 

|5| AT, 

y = ^^as,iXs,i (1) 

s=l 1 = 1 

where all the operations are performed in a Galois field of size q, GF{q). The ^-th original source packet of the s-th session 
is denoted as Xg^i and Ug^i is the corresponding coding coefficient. It should be noted that not all of the sessions necessarily 
participate in a particular inter-session network coded packet. When some of the sessions are not included in the combination, 
the corresponding coding coefficients for these sessions are zero. From that perspective, intra-session network coded packets 
can be viewed as a special case of inter-session network coding where packets only from one session participate in the coding 
operations. 

Depending on the available set of packets at the parent peers, every peer node may request from its parent peers intra-session 
network coded packets of its session, as well as inter-session network coded packets, i.e., packets that are combinations of 
different sessions. These combinations do not necessarily involve packets from the session requested by the peer 

'Scenarios where a user demands more than one session are not studied in this paper 
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All the possible packet types that can be generated in the network constitute the set T- Every element t of T represents 
a particular combination of sessions. Hence, in a network with |iS| concurrent sessions, the number of different packet types 
is 2''^ I — 1. Since intra-session network coding is a special case of inter-session network coding, intra-session network coded 
packets are also included in the set T- We denote as 7t the set of packet types that can be combined to generate packets of 
type t. The sessions that participate in a particular combination of packets t form the set St- We will refer to the sessions 
in the set St as the component sessions of flow of type t. We also define the sets and 7t,s. The set T" is a subset of 
T and contains the packet types that have session s as a component session, i.e., = {t £ T : s r\ St 0}. The set Tt.s 
includes all the packet types t' € Tt that can be used to generate packets of type t and have s as a component session, i.e., 
7t,s — {f (z Tt ■ s n Sf 0}. For example, if we consider packets that are combinations of the session Si and Sj, i.e., 
t = SiSj, then 7t = {si, Sj, SiSj}, St = {si,Sj}, 7t,s, = {si,SiSj} and 7t,sj = {sj,SiSj}. 

Every peer node, upon receiving a sufficient number of network coded packets, decodes the received packets in order to 
obtain the packets of the requested session. The decoding of a particular session is typically performed by means of Gaussian 
eUmination when a full rank system of packets is received. Note that, since the coding coefficients are drawn randomly 
according to a uniform distribution from the GF((7), a header of length J^ses ^^sio) bits is appended to the network coded 
packets. This header identifies all the coding operations performed on the packets while they travel through the network; it 
renders the decoding process feasible, since the encoding structure becomes implicit. 

In general, the application of inter-session network coding is not trivial. In fact, random mixing of all the available sessions 
is not always efficient, as it may cause an unacceptable increase of decoding delay for a specific source data. This is due to the 
fact that users need to receive enough innovative packets in order to decode all the encoded sessions along with the session of 
their interest. The term "innovative" refers to packets that bring novel information with respect to the packets that have been 
previously received by the peer These packets are linearly independent from the packets that are already stored in the peer's 
buffer In order to alleviate the shortcomings of the random mixing of all the sessions, an efficient rate allocation algorithm 
is essential. The goal of such rate allocation algorithm is to determine the sessions that should be combined and the rate that 
should be allocated to each combination in order to minimize the average decoding delay. This decoding delay depends on 
the innovative packet rates that the user receives for each of the session combinations that are available in the network. Let us 
denote as the input rate allocation vector associated with peer i, where each element of the vector represents the innovative 
rate of packets of type t & T that the peer i receives. The optimization problem that we aim to solve can be formally written 
as 



where Ai{ri,gi) is the expected delay needed for peer i to decode the source data gi £ S that peer i requests, and \Af\ is 
the number of peers in the network. The minimization problem in Eq. (|2]i seeks for the optimal rate allocation among all the 
possible session combinations in the network and determines the optimal network coding decisions in the network peers. It 
however necessitates the existence of a central control server that monitors the network conditions and has full knowledge of the 
network statistics. This is not practical in networks characterized by dynamics such as bandwidth variations, varying channel 
conditions, peer arrivals/departures at random time instances, etc. Therefore, we propose to optimize the decoding delay locally 
in a small neighborhood that comprises the peer itself and its parent and children peers. The rate optimization is performed 
locally with only a partial knowledge of the network statistics and the required communication overhead is minimal. Due to 
the distributed nature of the problem, global optimality can however not be guaranteed anymore, but the solution proposed 
below proves to be effective and adapted to realistic settings. 

B. Communication protocol 

The distributed delay optimization solution requires some exchange of information between the peers. We propose the 
following communication protocol. Let us consider the peer i and its local neighborhood that consists of the set of parent peers 
Ai and a set of children peers Vi as depicted in Fig. |2] We assume that the peer i is aware of the local network statistics, 
i.e., its total upload and download bandwidth ( and Cf, respectively), the channel capacity and the loss rates of the input 
and output links (bki, Tfki, Vfc e Ai and bij, iTij, Mj £ Vi, respectively). Every child peer j communicates to the peer i the 
identity gj of the session it wants to receive and the total download bandwidth C^. 

Whenever the peer i wants to optimize the requested packet flow rates, it requests the peers in its neighborhood to provide 
all the necessary information about the local status of the system. Specifically, every parent k G Ai sends to the peer i a 
vector Rk with the values of the input innovative flow rates for every packet type t € T. Every element of this vector 
represents the total input innovative flow rate of packets of type t available at the parent peer k at the time instant when the 
peer i performs the optimization of the rate allocation. In more details, i?^ is given as 




(2) 




(3) 
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Fig. 2. Communication protocol. Tlie neighborhood of the peer consists of the parent peers (in green) and the children peers (in blue). The green arrows 
(solid) indicate the information communicated to the peer i by its parent peers. The blue arrows (dashed) represent the information received by the peer from 
its children peers. 



where r^j, is the innovative rate of packets of type t received by node k from its parent node 7i0 Along with this value, every 
parent peer communicates to the peer i the maximum available upload bandwidth cjlj on the link (fc, i). This bandwidth is 
computed by subtracting from the total upload bandwidth the bandwidth that is required to serve all the children peers of peer 
k (except for the peer i). In other words, the available bandwidth is written as 

veT>k\i 

where we consider that the bandwidth used by the transmitting peer corresponds to the innovative rate observed by the receiving 
peer, augmented by the rate of the lost data. Similarly, every child peer j G Vi forwards to the peer i a vector Rj\i, with 
Vt £ T, representing the total innovative input flow rate of packets of type t that the peer j receives from its parent 
peers, except for the peer i 

Rh= E (5) 

Finally, every child peer sends to the peer i the value of the maximum available download bandwidth cfj, which is computed 
similarly to the available upload bandwidth in Eq. (|4]i as 

Practically, cfj is the remaining available download bandwidth of the child peer j, after subtracting the bandwidth that is 
needed to receive the requested innovative packets from all its parent peers excluding the peer i. The communication protocol 
is illustrated in Fig. |2] 



C. Distributed delay minimization problem 

We are now able to formulate the distributed rate allocation problem that is solved independently in every network peer It 
consists in determining the optimal innovative rates that the peer requests from its parents so that the average expected delay of 
the peer and its children peers is minimized. The reason for considering the children peers in the rate allocation optimization 
performed by every network peer is to avoid selfish behaviors of the peer nodes. It is obvious that if the peer performed the 
rate allocation taking into account only its own delay, it would preferably allocate all its resources to intra-session network 
coded flows, as there would be no incentives for the peer to request combined packets. In that case, the network peers would 
be unable to benefit from inter-session network coding. On the contrary, by including the delay of the children peers in the 
optimization objective, we provide incentives for peers to combine packets of different sources in order to serve as many 
peers as possible without major penalty on their own utility. By encouraging peers' collaboration we reach more socially fair 
solutions. 

Let us denote as r = r*^), Vfc G Ai, Vj G Vt £ T, the vector of innovative packet flow rates, where r|,j represents 
the innovative rate of packets of type t received by the peer i from its parent fc, while r*^ is the innovative rate of packets 
of type t received by the child peer j from the peer i. Formally, the distributed version of the delay optimization problem of 
Eq. (|2|i in the i-th peer node is stated as 

argininAj(r) s.t. r G 7?.™„ (7) 

r 

-Here we have assumed that two packets that anive from two different links are innovative with respect to each other with high probability. This holds in 
general in networks with high path diversity, which is the case considered in this work. 
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The search space TZmin is defined by a set of Hnear inequaHty constraints, which determine the set of feasible values of the 
innovative packet flow rates on the input and output links of the peer i 



ter 



< 2^ rl < 6fe,(l - ttm), Vfc e (8) 



< ^ 4 < b^Jil - 7T,j), Vj e V, (9) 
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The constraints appear in pairs and refer to the input and the output links of the peer i respectively. Eqs. dHJ and (|9]l are the 
link capacity constraints, which state that the sum of innovative packet rates for all packet types received on a link cannot 
exceed the hnk capacity. The second pair of constraints in Eqs. ( fTOb and (fTTI) asserts that the innovative flow after compensating 
for the losses cannot be larger than the download (upload) bandwidth of the node. Eqs. (fTSl) and ( fTSI ) reflect the fact that the 
peer cannot receive more than the maximum available upload bandwidth of the parent peers nor send more than the maximum 
available download bandwidth of the children peers. Eqs. (fT4b and ( fTSb give upper bounds to the innovative packet flow rates 
with the available innovative packet rates at parent nodes. Finally, the last pair of constraints in Eqs. (fTSI l and (fTTt limits the 
innovative packet rate by the available innovative rate provided by the sources, i.e., the peer cannot receive innovative packets 
faster than they are injected in the network by the sources. 

The average decoding delay of the peer i and its children peers Aj(r) is given by the following expression 

A.(r) = —i— (A,(r,gO+ E ^ji^^9,)) (18) 

The expected delay Ai{r, gi) experienced by the peer i for receiving and decoding a block of packets of the requested session 
gi depends on the average number of packets that the peer i needs to collect for decoding. The latter is a function of the types 
and the innovative rates of the packets that arrive at the peer 

The optimization problem stated in Eq. (|7]i is complex and in general non-convex. In order to solve it, we make the following 
simplifying assumptions. We assume that the time is slotted and that only one packet can be transmitted per time slot. We 
approximate the duration of the time slot by = Thus, we can estimate the average decoding delay as the product of 
the average time di required to receive one packet and the average number of packets E[l] that the peer receives before it is 
able to decode 

Adr,g,)^d,E[l] (19) 

The solution of Eq. ^ then requires the computation of the average number of packets E[l] that the peer and its children 
peers need to receive in order to decode their data of interest. In the following section, we will present an efficient method for 
computing E[l] that permits to transform the initial problem into a set of convex subproblems and to obtain a solution with 
low complexity. 

IV. Decoding delay analysis with equivalent flows 

The estimation of the average decoding delay as described in Eq. (fT9] l requires the computation of the expected number of 
packets E[l] that the peer has to receive for decoding one block of packets of the session of interest. The exact computation 
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of E[l] involves considering all the possible events that lead to a decodable set of I packets ||2TI . This is clearly non-trivial to 
compute. In this section, we introduce the notion of equivalent flows in order to approximate the decoding delay with simple 
functions that can be computed efficiently. 

Let us assume that the session s is the session of interest. There are several possibilities to decode the packets of this source 
data. Session s can be decoded from intra-session network coded packets when Ng such innovative packets are available. 
Otherwise, it can be decoded from a set of intra- and inter-session network coded packets of types t' e 7t, for any t e 7"^, as 
long as this set of packets forms a full rank system. In particular, when the decoding is performed from a session combination 
t, one needs A^^' innovative packets for each component session s' G 5*0 Note that these Ng' packets can be of any type 
t' e 7t,s'- Note also that any novel inter-session network coded packet of type t' E Tt contains novel information for all the 
component sessions s' G Sf- In other words, any novel inter-session network coded packet of type t' £ Tt can increase the 
rank of any of its component sessions. This property stems from the definition of innovation and is also guaranteed by the 
constraints ( fT4] i - ( fTTl i. 

The above observations bring us to the core idea behind the notion of equivalent flows. We can see that, when session s is 
decoded from the session combination t, we can treat every component session s' S iSt of t as a separate session for which 
we need to collect Ns' innovative packets of any type t' e Tt,s'- That means that the flow of packets of type t' can be split 
among its component sessions. The rate at which innovative packets are collected for the component session s' is equal to 
the sum of the contributions of each packet flow t' that has s' as a component session. The only difference between the case 
of decoding session s from intra-session network coded packets and the case of decoding the same session from the session 
combination t is that, in the latter case, the session s can only be decoded when a sufficient number of innovative packets is 
available for all the component sessions s' e St- We now propose a definition for the equivalent flows. 

Definition 1. Given a session combination t E T, we define an equivalent flow for every component session s E St of t as a 
virtual flow of packets with innovative rate equal to the sum of the contributions of innovative rates from every flow of type 

t' e Tt,s. 

We will henceforth refer to the rate of an equivalent flow as equivalent rate. Note that Definition [1] is general and applies 
also to types t that correspond to intra-session network coded packets. In this case, the equivalent flow coincides with the 
actual innovative flow of intra-session network coded packets. When t is a combination of two or more sessions, the innovative 
rate of the equivalent flow for every component session s E St is higher or equal to the actual innovative rate of the flow of 
intra-session network coded packets of this same component session. This increment in the equivalent rate comes from the 
contribution of the inter-session network coded packet flows that have the session s as a component session. 

Mathematically, the equivalent innovative rate vf ^ for the component session s in the session combination t received at the 
peer i can be represented as 

E <s E ^k^, E St (20) 

t'€Tt,s keAi 



where 7*^ E [0, 1] and J2s'es , X* s' — 1' ^ Tt- The coefficient 7*3 indicates the contribution of the innovative flow of 
type t' to the rate vj ^ at which innovative packets are collected for the component session s when session combination t is 
considered for decoding. 

Equipped with the definition of equivalent flows, we can now calculate the coefficients 7* ^ and approximate the decoding 
delay at peer i. Let us denote as 

V r* 

ft — fjd 

the probability of receiving an innovative packet of type t at peer i, where r^^ is the innovative flow rate of packets of type t 
that the peer i receives from its parent k. In a similar way, for a given session combination t E T and for every component 
session s E St, we define the probability 

t ^i.s J2t'eTt.s^hs'^keA^''^ki sr-^ t' t' /o-ix 

which represents the probability of receiving an innovative packet for the component session s at the peer i assuming that 
decoding is performed from the session combination t. The probability Pi sH) to receive the A^s-th innovative packet for the 
component session s of the session combination t upon receiving exactly I packets at peer i is given by the negative binomial 
distribution 

pul)-(^'^^i<llf^i^-cy-''' (23) 



'We say that a session s is decoded from the session combination t 6 T", when packets of types t' G 7t are used to decode the data of session s. 
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Thus, the average time needed for receiving Ng innovative packets for the component session s at peer i is 

°° N 

= d,E[l] = J2 ^PU^) = d~ (24) 

where E[l] stands for the average number of packets that the peer has to receive in order to collect Ns innovative packets 
for the component session s. It is given by the mean of the negative binomial distribution in Eq. (l23l l and, in our case, it is 
simply the ratio of the size of the block of source packets Ns and the probability ^ of receiving an innovative packet. Note 
that, in Eqs. ( l23T l and (|24] | we have assumed that the innovative rate remains constant with time. In practice, as the number of 
innovative packets in the peer's buffer increases, the probability of receiving a non-innovative packet also increases. However, 
in networks with high path diversity and for large Galois field sizes the probability of generating two identical or linearly 
dependent packets is negligible ll22l . This permits to make the assumption of a constant innovative rate. 

In order to determine the values of the 7* ^ coefficients for every session combination t we need to look at the problem 
from the point of view of the decoder The decoding of the session s from a session combination t is feasible as soon as Ns' 
innovative packets of any type t' e 7t,s' are available at the decoder for every component session s' £ St- This implies that 
the inter-session network coded flows are split among their component sessions in such a way that the delays for collecting 
the necessary number of innovative packets for every component session are as balanced as possible. That means that the 
equivalent rates, as seen by the decoder, are such that the maximum of the delays among the component sessions is minimized. 

We can now formulate the min max optimization problem that permits to determine the coefficients 7* ^ and subsequently 
the equivalent rates. The objective is to calculate the coefficients 7^ = {7|s} that minimize the maximum average delay A* ^ 
among the component sessions s ^ St. Formally, this optimization problem is written as 



At / A • .1 

mm max A, „(7i) — mmmaxdi— r — ; — - 

■y. sest sest qlsili) 

s.t. J2 76 <l'7t' e[0,l], WeTt 



(25) 



Once we have computed the equivalent rates, we can estimate the average decoding delay A* experienced by the peer i 
for decoding a block of packets of session s from the session combination t. Assuming that 7* is the optimal solution of 
the optimization problem in Eq. (l25l l, the decoding delay A* is simply the maximum of delays A* ^ over all the component 
sessions s ^ St- Indeed, in order to decode session s the peer needs to wait until the necessary number of packets is available 
for every component session. Thus, we have 

A* = inaxAt(7;) (26) 

Note that the vector of coefficients 7* is different for different session combinations t. 

To complete our analysis, we need to determine the average decoding delay observed at the peer i for decoding its session 
of interest s. Eq. (|26] | gives the average decoding delay under the assumption that the peer decodes from the specific session 
combination t that has s as a component session. In general, there may be multiple session combinations t' such that sDSf ^ 
yielding thus several possibilities for decoding. However, to simplify our analysis, we will assume that, for a given set of 
innovative packet flow rates on the peer's input links, the peer decodes the data of interest from the session combination that 
corresponds to the minimum average decoding delay A* 

Ai(r, s) = min A' = min max A* 5/(7*) (27) 

To summarize, in order to compute the approximate decoding delay for decoding packets of session s we have the following 
steps: 

1) we compute the equivalent flows by solving the min max problem in Eq. (l25T l for every session combination t e T*, 

2) using the equivalent flows computed in step 1, we calculate the approximate decoding delay for every session combination 
t from Eqs. ^ and (|26] |, 

3) finally, we approximate the delay with the minimum among the delays computed in step 2 (Eq. dZTli). 

Finally, it should be noted that we have considered the worst case scenario where all the component sessions involved in 
a session combination have to be decoded along with the requested session. This is due to the random encoding strategy 
deployed in our scheme. Other encoding strategies could be devised to avoid decoding all the sessions |23]. However, these 
strategies require expensive control and diminish the advantages of randomized network coding. The design of such encoding 
strategies is not trivial and is out of the scope of this paper. 

Example 1. We now illustrate the computation of equivalent flows and the estimation of the decoding delay with a numerical 
example. We assume that three sources, namely si, S2 and S3, are transmitted into the network and that the peer i requests the 
session si. The block sizes for the three sessions are Ns^ — Ns2 — Ns^ = 10 packets. We choose two sets of probabilities 
of receiving an innovative packet at the peer i for all the possible packet types t ^ T- These two sets of probabilities, shown 
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TABLE I 

Probabilities of receiving an innovative packet at the peer i for all the possible packet types. 
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TABLE II 

Probabilities associated with equivalent rates and the average number of packets needed for decoding. 
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0.2611 
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in Table U represent two different rate allocations at the peer i and correspond to two different instances of the decoding 
problem, namely Problem A and Problem B. Given these probabilities, we want to estimate the decoding delay at the peer 
i for decoding one block of packets from its source of interest. The session si can be decoded from any of the session 
combinations t — si,t — siS2,t — S1S3 ot t = S1S2S3. 

Table illustrates the results obtained by following the three steps summarized above. In particular, in Table |ll] we present 
the probabilities q* ^ that correspond to the equivalent flows for all possible session combinations that have session si as a 
component session. We also present the average number of packets that have to be received by the peer in order to decode 
its session of interest from the session combination t. Recall that the decoding delay is proportional to the average number of 
necessary packets for decoding (see Eq. (l24l). Finally, for comparison, we compute the average number of packets Ei with 
the exact method provided in [2 1 1 . 

From the results presented in Table HI] we can see that when the session si is decoded from intra-session network coded 
packets, the equivalent rate is equal to the actual innovative rate of intra-session network coded packets of type t = si, i.e., 
li^si — Pi'^ ■ On the other hand, when an inter-session combination is considered for decoding, the equivalent rates of the 
component sessions are higher than the innovative rates of the intra-session network coded flows of the component sessions. 
For example, when the session combination i = S1S3 is considered for decoding, we have > p^ ^ and ^^^,3^ > pl^- The 
increment in the rate comes from the splitting of the combined flow of type t = sis^ among its component sessions si and S3. 
Further, we can observe that in both Problems A and B the minimum number of packets required for decoding corresponds 
to the session combination t = 515253- We can also see that this number, calculated using the approach of equivalent flows, 
is very close to the actual average number of packets computed with the exact method liZTl . Another observation that we 
can make is that the performance in terms of decoding delay, for a given session combination, is driven by the component 
session that requires the most time in order to collect all the necessary innovative packets. Let us consider again the session 
combination t = S1S3. The equivalent rates for the component session S3 are almost the same in both Problems A and B, 
however, the equivalent rate for the component session si in Problem B is approximately half of the corresponding rate in 
Problem A and also half of the equivalent rate for the component session S3. Thus, in Problem B the peer needs to collect 
approximately two times more packets than in Problem A, in order to decode from the session combination t = S1S3. 

Finally, Fig. |3] illustrates the contributions of each packet flow to the equivalent rates of the component sessions si, S2 and 
S3, when session combination t = S1S2S3 is considered for decoding. Every color corresponds to a specific packet flow type, 
as shown in the figure. The bars represent the equivalent rates. The height of each bar is proportional to the magnitude of 
the corresponding equivalent rate. The height of a sub-bar of a certain color is proportional to the contribution of the flow 
denoted with the same color. We can see that, in Problem A, the packet flows are split among their component sessions in 
such a way that the equivalent rates for all component sessions are equal. On the contrary, in Problem B, we see that all the 
flows that have session si as a component session contribute only to the equivalent rate that corresponds to si. Furthermore, 
this equivalent rate is lower that the ones that correspond to the component sessions S2 and S3. 

V. Distributed rate allocation 

In this section we present the distributed rate allocation algorithm that solves the problem stated in Eq. ^ with the help of 
equivalent flow representations. According to our proposed solution, every peer solves a rate allocation optimization problem 
in two steps. In every optimization round, the peer first finds the optimal rate allocation that improves the average decoding 
delay of itself and its direct children. In order to find the optimal rate allocation, the initial problem is decomposed into several 
convex subproblems based on the equivalent flows representation described in Section |IV] Second, the peer maximizes the 
total throughput in terms of innovative packet rate while preserving the optimal rates obtained from the delay minimization 
step. This second step compensates for the partly myopic behavior of the peers and boosts the performance of the data delivery 
system as each peer can transmit packets that are potentially useful for other peers different than its direct children. 
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Fig. 3. Illustration of the contribution of each packet flow to the equivalent rates of the component sessions si, S2 and S3, when session combination 
t = S1S2S3 is considered for decoding. 



A. Decoding delay minimization 

The first step of our algorithm consists in finding the rates that minimize the decoding delay of a peer and its direct children. 
In order to determine these rates, the peer first obtains all the necessary information from its neighborhood following the 
communication protocol described in Section UlI-BI It then solves the rate allocation problem independently of the other peers 
and without any centralized control. 

The decoding delay minimization problem is stated in Eq. (|7]l. Recall that we have made a simplifying assumption that, for 
a given rate allocation, the peer i and its children peers j, (j e decode the requested data from the session combinations 
that correspond to the minimum decoding delay (see Eq. (l27t). Hence, the original problem can be decomposed into a set 
of convex subproblems. Every subproblem corresponds to finding the optimal rate allocation vector r = (r^,^ , r*^ ) , Vfc £ 
Ai,Vj e VijVt G T, that yields the minimum average decoding delay Ai(r) for a specific tuple of session combinations 
{U,{tjii G T^i}) e T^' X J] T^J. Combining Eqs. (|7]i, ([TSll, (|26ll and (|21l, the subproblem of finding the optimal rate 

allocation for a given tuple (t^, {tj, j G G 7"^' x 11 T^' , can be written as 



arffmin \ (di max —. — (- > dj max —. — - — ) 

(28) 

\ " _ t' 



2^ 7L' <l, Vt'e7I„, Vne{zUP,} 



l.S 



From Eq. (|22] |. q\ ,, = J2t'ert s ^i.sP\ ■ replace the product 7* ^p* with the variable x G [0, 1] and write q, 

J2t'£Tt s- minimization problem in Eq. ( |28] | becomes 

are; mm - — ; a, max —i h > di max — r 

s.t. r G7^™„ (29) 
^ x^;,, < pi, W e Vn e {^ U V,} 

where = — ^^^f — ^ and p* = ^ , j G C;. The problem stated in Eq. (|29] | is convex. This can be shown using 
the following arguments. The function t ^ is convex since it is a composition of the convex function - with the affine 

expression ql^^{xi). The pointwise maximum is also a convex function. Thus, the objective function in Eq. (|28] | is convex 
since it is a nonnegative weighted sum of convex functions [24,1 . It can be solved using the CVX Matlab-based package [251 
for example. 

The solution to the initial rate allocation problem stated in Eq. (|7]i can be obtained by solving the subproblems of Eq. (|29] l for 
all the feasible tuples (ti, {tj,j e 2?^}) e T^' x J| T^j . The results of these subproblems are then combined and the solution 
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(rate allocation vector) that yields the minimum delay is chosen. This solution also constitutes the solution to the original 
problem in Eq. O. The number of convex subproblems to be solved depends on the cardinality of the set T^' x Y[ 

that grows exponentially with the number of sources available in the network and the number of the peer's direct children. In 
practice, however, the number of sources is typically small, and the peers have a limited upload bandwidth, which allows only 
a few children nodes to be connected simultaneously to the same peer. Therefore, the number of convex subproblems to be 
solved by each node is typically small. 



B. Maximization of the total innovative input rate 

The solution of the minimization problem in Eq. O guarantees the transmission of data sources that are requested by the 
peer i and its children peers at optimal rates, as long as these data sources are available at the parent nodes of the examined 
peer. However, a peer is not aware of the data requested by peers that are two or more hops away. This means that, in certain 
cases, some of the sessions that are available in the network and may be potentially useful for other peers beyond the peer's i 
neighborhood, are never requested by the peer i. Thus, these sessions can never be forwarded when requested by other peers, 
which eventually penalizes the performance of the downstream nodes. This drawback of the distributed scheme is a result of 
having a limited network horizon with only local information in solving the rate allocation. In order to reduce the effect of 
this shortcoming, we propose to solve a simple throughput maximization problem. This maximization problem is solved in 
every optimization round immediately after the optimal rates have been determined as presented in Section IV-AI Specifically, 
we aim at maximizing the total innovative packet flow rate for all the packet types such that the flow values are larger or equal 
to the optimal flow rates computed from Eq. (|7). Practically, this means that, whenever there exists some unused bandwidth, 
it is allocated to packet flows that are not explicitly requested by a peer or its children peers, but that can be potentially useful 
for other nodes. The maximization problem can be formally written as 

argmax ^ ^'rl^ 

keA, ter (30) 



where r^^^ = argmin Ai{r) and the inequality sign between two vectors denotes the inequality relationship between vector 

r 

elements at the same positions. The search space TZmax is given by the following set of linear inequality constraints 



< J2 ^ - ""ki), Vfc e A, yter (31) 

EEn^.^^t (32) 

5]rL <c)^,(l-7r,0, Vfce A (33) 
ter 

^k^< E ^k, Vt e r, V.S e St, Vk e A (34) 

Y'^k^< Us, V< e r, Vs e St (35) 

t'STt,, keAi 

This set of constraints has a similar interpretation as the one that defines the search space for the delay minimization problem 
(Eqs. (O to ( [Tt] )) except for the fact that only the constraints of the input links are considered. The optimization problem stated 
in Eq. ( l30t is a linear program and can be solved using any of the standard optimization algorithms Ii24l . 

Note that, since r^^ is the rate at which innovative packets of type t arrive at peer node i from its parent peer fc, the actual 
rate /^^ that the peer i has to request from its parent should be augmented by the average packet loss rate that is observed on 
the link 

fL = (36) 

-1- — TT/ci 

The communication protocol and the distributed rate allocation algorithm are summarized in Algorithm [T] The algorithm 
runs periodically in every peer This allows to adapt the rate allocation to possible changes that may occur in the network. In 
practice, a peer optimizes its input rates only when all its parent peers have also performed the optimization. The optimization 
stops when the utility of the peer does not change for a certain number of optimization rounds and all its parents have stopped 
optimizing or if the peer has reached a maximum number of optimization rounds. 
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Algorithm 1 Distributed Rate Allocation Algorithm. 
1: Initialization 

2: Set the current optimization round li — {li — oo for source peers). 

3: Define the maximum number of optimization rounds Imax, the minimum number of optimization rounds Imin and the 
number of optimization rounds Is- 

4: while It < Irnax dO 

5: Request the values Ik, V/c G Ai, from the parent peers. 
6: a Ik > k yk G Ai then 

7: Request the values of and c^!^, \/k e Ai, from the parent peers and the values Rj\i, cfj,gj and Cj, Mj E Vi, 
from the children peers. 

8: Solve the delay minimization subproblem (Eq. ( |29] l) V(ti, tj,jeVi ) G 7"^* x H ■ Combine the results and determine 

the optimal rate allocation vector r = {rki,rij),\/k G Ai, Vj G Vi. 
9: Solve the throughput maximization problem (Eq. (|30] |) and update the rates r^^ , Vfc G 

10: Compute the actual rates fli^Vk G yli to be requested from the parent peers, as /^^ = jz^^ 

11: if Ik —— oo Vfc G Ai and > Imin and Ai(r, g^) has not changed for Ig rounds then 

12: Set li — oo 

13: else 

14: Set li ^ li + 1. 

15: end if 

16: else 

17: Go to step 4. 

18: end if 

19: end while 



VI. Performance evaluation 

We analyze the performance of our proposed distributed inter-session rate allocation algorithm. We consider the transmission 
of data from multiple concurrent sessions and evaluate the performance of the proposed scheme in terms of the average decoding 
delay. The decoding delay is measured as the time needed for a network peer to collect and decode one block of packets from 
the source of interest. 

First, we provide an in-depth study of the behaviour of our rate allocation scheme in a small size toy network. We then present 
the result of applying the proposed method to larger topologies. We compare the performance of our scheme, henceforth denoted 
as "InterNC" (Inter-session Network Coding) to a baseline intra-session network coding rate allocation scheme "IntraNC" (Intra- 
session Network Coding). The latter is a modification of the proposed method except for the fact that the coding across different 
sessions in the network nodes is not allowed. For the sake of completeness, we also provide a comparison of the decoding delay 
and the optimal rate allocation with a centralized algorithm that solves the rate allocation problem as given in Eq. (|2). Note 
that the centralized scheme assumes full knowledge of the network statistics, and has a complexity that grows exponentially 
with the number of nodes in the network, so that it does not represent a viable solution in large networks. 

A. Toy network 

We first evaluate the performance of the proposed distributed inter-session rate allocation algorithm for the network depicted 
in Fig. Ufa). The network consists of 3 source nodes and 9 peer nodes, which subscribe to different sources. The packet loss 
rate is set to 5% on all links. The bandwidth of the links that originate from the sources, as well as of the link connecting 
nodes and ng is set to 30 packets/sec. The bandwidth of the links that originate from nodes nr,ns and ng is set to 60 
packets/sec. The block size for all 3 sources is 10 packets. 

Fig. HJb) presents the evolution of the average delay of the network clients with respect to the bandwidth of the links 
connecting nodes 714,717 and ne,ng for all the schemes under comparison. We can observe that, even for low link rates, the 
proposed distributed InterNC rate allocation scheme performs better than the distributed IntraNC scheme. The gains come 
from the fact that the nodes can combine packets from different sessions on bottleneck links, whereas in intra-session network 
coding the performance is limited by the presence of low rate links that cannot serve all the clients at the same time. As 
the link rates increase, higher gains in terms of delay can be noticed for our proposed InterNC scheme, as more packets are 
combined across different sessions. On the contrary, the IntraNC schemes fail to deal efficiently with the bottleneck created on 
the link between the nodes and 7is and the slight improvement of the average decoding delay comes only from the increase 
of the rate at which packets are supplied to node nu. 

Finally, we can notice that the distributed rate allocation schemes, both the proposed InterNC scheme and the baseline 
IntraNC network coding scheme, manage to reach the performance of their centralized counterpart. This essentially means that 
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Fig. 4. (a) Toy P2P network topology where 3 data sources ai'e concurrently transmitted to the peer nodes. The source label next to each peer indicates the 
source data that this node wants to receive, (b) Average decoding delay for the toy network topology depicted in Fig. |4ja) versus the bandwidth of the links 
connecting nodes n4 , and ne , ng . 



for this specific network topology the limited knowledge of the local network statistics that is available to the distributed rate 
allocation algorithms is sufficient to achieve the global optimal rate allocation solution that can be attained by the centralized 
schemes. However, we expect that in generic topologies the performance of the distributed rate allocation algorithms will be 
inferior to that of the centralized ones, as the myopic optimization performed by the peers does not always detect all the 
opportunities for inter-session packet combinations. 

Our conclusions regarding the average decoding delay can be further supported by examining the innovative rate that is 
achieved by the schemes under comparison. Figs. |5] and |6] illustrate the normalized total innovative input packet rate of nodes 
riT, ns and rig for the distributed and centraUzed algorithms, respectively. The normalization is done with respect to the total 
input bandwidth of the peer In the figure, Sj denotes a flow of intra-session network coded packets of session Sj, whereas 
SiSj represents the combined flow of inter-session network coded packets from sessions Si and sj. The flows that are zero in 
the whole range of link bandwidths are omitted from the figures. 

As we can notice from Figs. |5lb) and |3b), the link between nodes and rig has to be shared by the flows si and S3 
when only intra-session network coding is allowed, as this is the only path from where nodes riio and ni2 can receive their 
requested flows. Thus, when the bandwidth of the Unks between nodes ri4, n-j and ng, ng increases, the average decoding delay 
of nodes riio and ni2 cannot be improved as they receive intra-session network coded packets at constant rates regardless of 
the bandwidth variations. The only reason for the slight improvement of the average delay that we observe in Fig. 2Ib) is the 
additional supply of packets of session S2 to node nn from node nj, as can be seen by observing the rate curves in Figs. |3a) 
and|6ja). 

When inter-session network coding is allowed, the average performance of the network is enhanced mainly by the combination 
of flows si and S3 on the bottleneck link between nodes and uq. As we can see in Figs. |5lb) and|6jb), the node ng allocates 
part of the input bandwidth to the combined flow S1S3 whereas the rest is allocated to the intra-session network coded flow 
S3. As the node ng starts to provide more intra-session network coded packets of flow S3 to node ni2 when the bandwidth 
increases, the percentage of rate for the combined flow on the bottleneck link increases and eventually the node ng requests 
only combined packets. At this point, both nodes nio and ni2 manage to receive their requested flows at the rate of the 
bottleneck link since they receive at the same rate the other component packets of the combined flow from nodes riy and ng 
respectively, and they are able to decode faster the session of their interest. Thus, we can see that the limitations imposed by 
the bottleneck link can be overcome by deploying inter-session network coding and utilizing the additional resources of the 
nodes for receiving packets that can help in decoding the combined sessions. 

It is worth noting that the rate allocation achieved by the distributed inter-session network coding algorithm is not identical 
to the one achieved by the centralized scheme for link bandwidth equal to 30 packets/sec, as can be seen in Figs. |3b), [Sjc), 
|6jb) and |6jc). This is attributed to the fact that the centralized algorithm has the full knowledge of the network topology. It 
can detect more opportunities for combining packets from different sessions, whereas the distributed scheme can only take 
advantage of the local network conditions. Finally, we can observe that for all schemes, the innovative rates and the average 
delay saturate as links' bandwidth reaches the value of 60 packets/sec. This is essentially the point where the system has 
reached the state where no other improvement can be achieved with either of the schemes. 

Note that, since nodes ^4, and ng receive all their packets directly from the sources, they do not affect the average 
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Fig. 5. Normalized total input innovative packet rate for nodes (a) n7, (b) ng, (c) ng versus the bandwidth of Hnks connecting nodes n4,n7 and ngjng 
for the topology depicted in Fig. |4}a). The schemes under comparison are the distributed InterNC and the distributed IntraNC rate allocation algorithms. 
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Fig. 6. Normalized total input innovative packet rate for nodes (a) n7, (b) ng, (c) ng versus the bandwidth of Hnks connecting nodes n4,n7 and n6,ng 
for the topology depicted in Fig. |4}a). The schemes under comparison are the centralized InterNC and the centralized IntraNC rate allocation algorithms. 

observed delay. The behavior of nodes nio, nn and ni2 depends also on the rates available at nodes ny, and ng as by the 
construction of the network they have sufficient download bandwidth in order to download all the packets that are available 
in the aforementioned nodes. 



B. Larger networks 

In this set of experiments, we evaluate the performance of the proposed scheme for the clustered network depicted in Fig. [T] 
This network consists of three server nodes and 30 client nodes. The clients are organized in 3 clusters of 9, 12 and 9 nodes 
respectively. Each cluster is an irregular directed network generated from a regular network by removing and shifting randomly 
some of the links ||26l . The pruning and shifting probabilities are set to 40% and 20% respectively. Every peer node is assigned 
one of the data sources. The selection of the sources is done uniformly at random. The clusters 1 and 3 are connected directly to 
the servers with links that have a capacity of 468 kbps each, whereas the cluster 2 is connected to the clusters 1 and 3 through 
links with a capacity that varies in the interval [117, 702] kbps. Moreover, the cluster 2 receives some packets directly from the 
sources through low speed links that have capacity of 117 kbps. Finally, the nodes within all the clusters are interconnected 
with high speed links of 1.6 Mbps. The packets size is fixed to 1500 bytes including the network coding header Again we 
consider that the block sizes for all data sources are equal to 10 packets. All the results in this section are averages of 10 
random realizations of the network. 

Fig. [HI a) illustrates the average decoding delay for the clustered network depicted in Fig. |7] with respect to the bandwidth 
of the links that connect cluster 2 to clusters 1 and 3. The schemes under comparison is the proposed distributed InterNC rate 
allocation algorithm and the baseline distributed IntraNC scheme. We can observe that, by allowing peers to combine data 
from different sessions, we can achieve lower decoding delay times than those that can be achieved with intra-session network 
coding only. As presented in Fig. |8lb) the gain is observed in cluster 2 that does not have sufficient resources to provide 
intra-session network coded packet to all the peers, contrarily to clusters 1 and 3 where all the peers are able to acquire all the 
packets directly from the sources. Thus, inter-session network coded packets are requested on the bottleneck links connecting 
cluster 2 to clusters 1 and 3 in order to serve more peers in the network, whereas the additional packets that are provided 
through the low capacity links that connect cluster 2 to the sources are used to decode faster the combined packets. 
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Fig. 8. Performance comparison of tlie InterNC and IntraNC algorithms regarding the average decoding delay as a function of the links' bandwidth for the 
examined cluster network: (a) average over all the nodes in the network and (b) average for each cluster of the network separately. 



VII. Conclusions 

We have proposed a novel distributed rate allocation algorithm for delivery of multiple concurrent sessions in wireline mesh 
networks. The algorithm is based on inter-session network coding. The network peers decide locally on the optimal coding 
decisions and rates for each combination of packets that they request from their parent peers. The decisions are based on 
the minimization of the average decoding delay of the peer and its children peers and require only a minimal communication 
overhead. We show that the initial non-convex rate allocation problem can be decomposed into a set of simpler convex problems 
with the help of the new equivalent flow representation. The final rate allocation can then be obtained by combining the results 
of each of the subproblems. The evaluation of the proposed algorithm demonstrates the benefits of utilizing inter-session 
network coding in terms of the decoding delays and efficient exploitation of network resources. The simulations also show 
that our distributed rate allocation algorithm performs closely to a centralized rate allocation solution in small networks. Our 
algorithm therefore offers a constructive solution for the design of effective inter-session network coding systems in dynamic 
networks. 
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